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Abstract 

In this paper we investigate the various aspects of noise and order in the 
micromaser system. In particular, we study the effect of adding fluctuations 
to the atom cavity transit time or to the atom-photon frequency detuning. By 
including such noise-producing mechanisms we study the probability and the 
joint probability for excited atoms to leave the cavity. The influence of such 
fluctuations on the phase structure of the micromaser as well as on the long- 
time atom correlation length is also discussed. We also derive the asymptotic 
form of micromaser observables. 
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1 Introduction 



Noise is usually considered as a limiting factor in the performance of a physical device 
(see e.g. Refs.0). There are, however, nonlinear dynamical systems where the presence 
of noise sources can induce completely new regimes that cannot be realized without noise. 
Recent studies have shown that noise in such systems can induce more ordered regimes, 
more regular structures, increase the degree of coherence, cause the amplification of weak 
signals and growth of their signal-to-noise ratio (see e.g. Refs. |2|, [§). In other words, 
noise can play a constructive role, enhancing the degree of order in a system. 

The micromaser is an example of such a nonlinear system. The micromaser system is 
an experimental realization of the idealized system of a two-level atom interacting with 
a second quantized single-mode of the electromagnetic field (for reviews and references 
see e.g. Refs. Q- J7|). In the micromaser a beam of two-level atoms is sent through a 
microcavity where each atom intersects with the photon field inside the cavity during a 
transit time r. After exit from the cavity the atoms are detected in either of its two states. 
It is assumed that subsequent atoms arrive at time intervals which are much longer then 
the atom-field interaction such that at most one atom at a time is inside the cavity, which 
is the operating condition for the one-atom maser. In such a system noise-controlled 
jumps between metastable states have been discussed in the literature ||. 

In the present paper we study the effect of including noise-producing mechanisms 
in the micromaser system like a velocity spread in the atomic beam or a spread in the 
atom-photon frequency detuning Aw. We show that under suitable conditions, such noise- 
producing mechanisms lead to more pronounced revivals in the probabilities V(+) and 
V(+,+), where V(+) is the probability that an atom is found in its excited state after 
interaction with the photon field inside the microcavity and V(+, +) is the joint proba- 
bility that two consecutive atoms are measured in their excited state after interaction. It 
is the purpose of the present paper to study the counterintuitive role that noise can play 
in the micromaser and extend the results in Refs. || § . 

The paper is organized as follows. For the convenience of the reader we recapitulate 
in Section [2] the theoretical framework of the micromaser system. In this Section we also 
discuss the general conditions to make revivals in the probability V(+) well separated 
in the atomic transit time. In Section |^ we determine asymptotic limits of the order 
parameter of the micromaser system and the probabilities V(+) and V(+, +). In Section 
^ the effect of noise on V(+) and V (+,+), the order parameter and the micromaser phase 
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diagram is discussed. In Section |5| the effect of fluctuations on the correlation length is 
discussed and numerical investigations are presented. A conclusion is given in Section []. 



2 The Dynamical System 



In the description of the dynamics of the one-atom micromaser we have to take losses 
of the photon field into account, i.e. the time evolution of the photon field is described 
by a master equation. The continuous time formulation of the micromaser system [10 is 
a suitable technical frame work for our purposes. Let a be the probability for a pump 
atom to be in its excited state. Assuming that the pump atoms are prepared in an 
incoherent mixture, i.e. the density matrix of the atoms is diagonal with diagonal matrix 
elements a and b such that a + b = 1, it is shown in |K| ||] that the vector p formed by 
the diagonal density matrix elements of the photon field obeys the differential equation 
dp/dt = —'-yLp. Here 7 is the damping rate of photons in the cavity and L = Lc — 
N(M - 1), where (L c ) nm = (n& + l)[n5 n>m - (n+ l)5 n+ i im ] + n b [(n + l)8 n>m - n5 n , m+ i] 
describes the damping of the cavity, n& is the number of thermal photons and N is the 
average number of atoms injected into the cavity during the cavity decay time. The 
matrix M is M = M(+) + M(-), where M(+) nm = bq n+ i5 n +i,m + a(l - q n +i)5 n ,m and 
M(—) nm = aq n 5 ntm+ i + b(l — q n )8 n ,m, has its origin in the Jaynes-Cummings (JC) model 
0, M. The quantity q n = q(x = n/N) is given by 



= ^^in 2 (V^+A2) (1) 

where A 2 = (Au) 2 /(Ag 2 N) is the scaled dimensionless detuning parameters of the mi- 
cromaser and g is the single photon Rabi frequency at zero detuning of the JC model. 
Eq. ([I]) is expressed in terms of the scaled dimensionless pump parameter 9 = gry/N. 
Noise mechanisms, like a spread in the transition time or a spread in the detuning, are 
now included in the analysis by simply averaging the matrix L with respect to 9 or A, i.e. 
averaging the quantity q n with respect to one of these parameters. A similar averaging 
procedure with regard to the parameters a (or b) or n& leads only to a trivial replacement 
of the corresponding parameters with their mean values. The stationary solution of the 
photon distribution where such noise effects are included is therefore derived in a standard 
and well known manner || [10, ||. The result is 
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nn b m + Na(q m ) 
m =i ^ + n b) m + Nb{q m ) 
where (•) in Eq. (||) denotes averaging, to be discussed below, with respect to 6 or A. p 
is a normalization constant. 

After the passage through the microcavity we make a selective measurement of the 
atoms. We imagine that one then only measure those atoms leaving the cavity with a 
definite value of 9 or A, i.e. in effect putting a sharp velocity filter or a filter sensitive to 
detuning at the atom output port of the cavity. The probability V(s) of finding an atom 
in a state s — ± after the interaction with the cavity photons, where + represents the 
excited atom state and — represents the atom ground state, can then be expressed in the 
following matrix form ||: 

V{s) = u T M(s)p , (3) 

such that V(+) + V{—) = 1. The elements of the vector p is given by the equilibrium 
distribution in Eq. (|2]). The quantity u is a vector with all entries equal to 1, u n = 1. 
Explicitly V(+) takes the form 

oo oo 

P(+) =a ^2 ~ qn+l ^ + b Pn+lQn+l ■ (4) 

n=0 n=0 

It is well known in the literature that this probability can exhibit quantum revivals (see 
e.g. Refs. Furthermore, the joint probability for observing two consecutive 

atoms in the states s\ and S2 is given by 



V{si,s 2 ) = u T S(s 2 ) S(si) p 

= u T M(s 2 ) S( Sl ) p , (5) 
where S(s) = (1 + L c /N)- l M(s). Explicitly one finds 

oo 

V(+,+) = a 2 (1-^+0(1 + L c /N)-^(l-q m+1 )p m 

n,m=0 



oo 

+ abJ2 { (1 - ?n+i) (1 + L c /N)-l 

Qm+l Pm+1 



n,m=0 
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+ q n (l + L c /N) n ^(l-q m+1 )p r , 



oo 



+ b 2 } q n (1 + L c /N) n ^q m+1 p m+1 . (6) 



n,m=0 



In Refs.|| [12| it is shown that the equilibrium distribution Eq. (0) can be re- written 
in a form which is rapidly convergent in the large N limit by making use of a Poisson 
summation technique |13[]. In terms of the scaled photon number variable x = n/N the 



stationary probability distribution can be written in the from 



where 

oo 

V(x) = V k (x) . (8) 

k=— oo 

The effective potential V(x) is expressed in terms of 

V k (x) = -[ dv \n[w(v)} cos(2n Nkv) , (9) 

where 



o 



n b x + a{q(x)) 

W[X) = — r —. — -, — tt . (10) 

V ; (1 + n b ) x + b (q(x)) V ; 

We stress that Eq. (|7|) is exact. In the large N limit Eq. (|7|) can be simplified by mak- 
ing use of a saddle-point approximation. Apart from calculable 1/N corrections we put 
V(x) = Vo(x). The saddle-points are then determined by Vq(x) = 0, where Vq(x) is the 
derivative of Vq(x) with respect to x. Hence, it is the nature of the global minima of Vq(x) 
which determine the probability distribution, apart from possible zeros of w(x). 

If the only global minimum of Vq(x) occurs at x = 0, which corresponds to the thermal 
phase of the micromaser, we can expand the effective potential Vo(x) around origin. The 
probability p n is then given by 

n b + a 8 2 



eff 

n b + b6 2 eff 



Pn = Po \ ; , z . rm 1 > (11) 



where 
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(q(x)) _ /sin 2 (#A)\ 



x^O X 



9i ff = limX±^- = {—± r J-) . (12) 



A 2 / 



The probability distribution as given by Eq. flTTD is normalizable provided 

9 2 eff (2a-l)<l . (13) 

If, on the other hand, there exists non-trivial saddle-points of Vq{x), which correspond 
to the possible maser phases of the micromaser, we write 

p( x ) = ^2pj( x ) > ( 14 ) 

3 

where Y2j denotes the sum over the local minima of Vq(x), i.e. Vq(x) > 0, and where 
pj(x) is p(x) for x close to the local minimum x = xj. The distribution pj(x) is given by 

p i(x ) = Jj= e -f vS^X*-*tf , (15) 
v 2irN 

where Tj is determined by the normalization condition for p(x), i.e 

e -NV (xj) 

Tj = e^^Vv^CI ' (16) 

and 

r /f f M _ ( 2q ~ l ? jq{x)) -xd(q(x))/dx 

V ° (X) - a + n b (2a-l) x* ' [L7) 

For the given parameters, the sum in Eq. (|16D is supposed to be taken over all saddle- 
points corresponding to a minimum of V (x), i.e. all saddle-points corresponding to 
Vq(x) > 0. If x = Xj does not correspond to a global minimum of the effective potential 
Vq(x), then Tj is exponentially small in the large iV limit. If x — Xj does correspond 
to one and only one global minimum, we can neglect all the terms in the sum in Tj but 
m — j, in which case Tj is reduced to Tj = \/VqJxJ). In the neighbourhood of such a 
global minimum the probability distribution in Eq. fll~5|) is therefore reduced to 



Pj{x) 



where the standard deviation o x of Eq. (|T8"D is 
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which is zero in the large N limit provided V^Xj) 7^ 0. The probability p(x) is therefore 
peaked around the global minima for large N. The minimum Xj of the potential Vq(x) is 
a solution of saddle-point equation 

x = (2a-l)(q(x)) < 1 . (20) 

This saddle-point equation is independent of the number n& of thermal photons in the 
cavity. 

The excitation probability V(+) can be re-written in a form where quantum revivals 
(see e.g. Ref. 0) become explicitly by making use of a Poisson summation technique 
p~3|j . When a = 1 we obtain 



1 OO -1 CO 

s 1 V - ^ 71 + 1 1 \ ^ 



2^^ n + l + iVA 2 2 

n=0 1/=— 00 



where 



1 Po TTna* C0S ( 2d VW+^) , (21) 



w (6) = N l°°dx p(x) - lylj*^ cos (2 Qyjx + 1/iV + A 2 ) , (22) 



and 



X r ^ P(X) x+l/N + A* V ^ m ^- d ' 2 ^ 2 } ■ (23) 

Here p(x) is the continuous version of p(n/N). We stress that Eq. ( PTj ) is exact. If a ^ 1 
we obtain 7- > (+) by the replacement P(+) — > 1 — a + (2a — 1)V(+), apart from 1/iV 
corrections. If the probability distribution p(x) is sufficiently peaked, i.e. x ^> a x , where 
x = x(8) denotes the average of x = n/N with respect to the stationary probability 
distribution p n , then the excitation probability is reduced to 
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1 ~~ 1 

P(+)Kl -2FTAJ + 2E^) . (24) 

u=0 



where 



x 

W u (9) Hi p(x = x v ) 



x + A 2 7rv / 2lA/v 



X COS 



9 x2 



2vrz/iV + A . 

2nuN 4 



7l 



v > 1 , (25) 



and 



H^) 2 - A! ' (26) 

According to Eq. fl26|) the i/:th revival of 7 ? (+) occurs in the region where 9 is close to 

NVx + A 2 . (27) 

This equation relates the width a x of the probability distribution p(x) into a measure 
for the width A9 U in the pump parameter of the v.th revival according to a probability 
distribution of the form p(9) ~ exp(— (9 — 9 l/ ) 2 /2(A9 u ) 2 ), i.e. 

A0„ = vrWV - ffa . (28) 

Two consecutive terms of the sum in Eq. (E3J) separate in time when their temporal 
separation 6 V+1 — 9 U is larger then A9 V+1 + A^, i.e. within one standard deviation, 
provided 

x + A 2 1 

^< 5 • (29) 

The revivals in P(+) cannot be resolved for values of v larger then the right-hand side of 
Eq. 

If the probability distribution p(x) is sufficiently peaked, i.e. x ^> a x , then the excita- 
tion probability V(+) approaches 

V(+) = a - i(2o - 1) (^ ( _ + X A2) - «*>(0)J w o - (2o - l)g(x) , (30) 
in the large N limit, where we make use of 
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w (0) w exp(-e 2 /2N) X - - cos(2flVx + A 2 ) . (31) 

x + A 2 



If the average in the saddle-point equation Eq. (j20|) is such that (q(x)} « ?(^), we see 



that 7- > (+) ~ a — x. If we average Eq. (|3"T|) with respect to the noise in the system, we 



obtain (V(+)) ~ a — x by making use of the saddle-point equation Eq. (|20D once more. 

The approximative result Eq. ([3(]) can actually be converted into an exact and more 
general relation between x and V(+). If the atoms are not measured with well-defined 
9 or A then we must in general average V(+) with respect to the corresponding proba- 
bility distribution. In Eq. (^) we then replace q n by (q n ) and denotes the corresponding 
excitation probability V(+) by (V(+)). From the equilibrium distribution Eq. (0) we 
then obtain the recurrence formula [ (1 + n&) n + Nb(q n ) } p n = p n _i[nb n + Na{q n ) ]. By 
summing this formula over n, we therefore derive that 

x = a+^-(V(+)) . (32) 

In general and in the presence of noise in 9 and/or A, x and V(+) have no direct relation 
and we must in general consider them to be independent variables. 



3 Absence of Fluctuations 



If the atoms enter the cavity without a spread in the transit time and without a 
spread in the atom-photon frequency detuning, the behavior of the order parameter and 
the excitation probability V(+) is well known in the literature (see e.g. Refs. @, || ||). As 
e.g. seen in Fig. [l], V(+) under such circumstances shows no clear evidence for resonant 
behavior of revivals. A more precise way to express the presence of revivals in the temporal 
variations of V(+) can e.g. obtained by performing a Fourier transformation of V(+) 
and considering the width of the corresponding spectrum. The corresponding Fourier 
spectrum of Fig. [I] is then broad. The appearance of clear revivals would correspond to a 
more peaked Fourier spectrum. 

In the large and 9 limit the order parameter approaches a constant. This constant, 
Xoo, is therefore determined by the global minimum of Vq(x) in the large 9 limit, in which 
case Vq(x) has a unique minimum. The micromaser system has then no maser-maser phase 
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transitions. As shown in an Appendix, this minimum is determined by dVo(x)/dx = 0, 
where 



+ f(y(x))-f(w(x)) , (33) 

and where f(z), y(x) and w(x) are as given in the Appendix. When e.g. a — 1, n& = 0.15 
and A = the global minimum of Vq(x) occurs at Xoo ~ 0.34, which corresponds to 
an asymptotic excitation probability V 00 ( J r) ~ 0.66 due to Eq. (^). When a = 1 and 
A = it follows from Eq. @ that for large 9 the joint probability V(+, +) approaches 
the asymptotic value 

7>oc(+,+)~(5Poc(+)-l)/4 , (34) 

i.e. with the physical parameters as above we find Poo(+, +) ~ 0.57 (see Fig. [I]). Eq. ( [34] ) 
follows from Eq. (P) by making use of the following properties of the unbounded operator 
L c 

oo oo 
J2( L c)nm = , ^ (L C ) nm = -1 , (35) 
n=0 m=0 

and performing a suitable large iV limit. 



dV {x) 
dx 



In 



4 Effects of Fluctuations 



We will now consider physical effects of noise in the micromaser system. We start by 
discussing the effects of adding fluctuations to the pump parameter 9 and then, in Section 
4~2|, we study the effects of adding fluctuations to the atom-photon detuning parameter 



4.1 Pump Parameter Fluctuations 

Suppose that the pump parameter 9 is described by a positive stochastic variable £ 
as described by the probability distribution P#(£) such that (£) = 9 and ((£ — #) 2 ) = oj. 
The averaged value of q'(x) with respect to Pe(£) is then given by 
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Figure 1: The probabilities V{+) (solid line) and V(+,+) (dotted line) as a 
function of 9 when the g m -terms in Eq. (||) are not averaged. The parameters are: 
a = 1, A = 0, nt, = 0.15 and N = 35. The curves show no clear evidence for 
the resonant behavior of revivals. In the large 6 and N limit the probability V{+) 
approaches 7 ? 00 (+) ~ 0.66 and the probability V(+,+) approaches 7 3 0O (+,+) ~ 
0.57, which are indicated by solid lines in the figure. 



(q(x))e = J™ dZPeiO^^ sin 2 (tV^TAl) . (36) 
Non-trivial saddle-points of the corresponding Vq(x) can be found by solving the equation 

l = (2a-l)h(6,x(6)) , (37) 

where 

h(0, x) = J™ di P e (0 sin 2 (^v^+A^ ) < a 2 e + 6 2 . (38) 

For small values of x, we can expand the effective potential V${x) around x = 0. A 
straightforward expansion Vq(x) then leads to 
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Vq(x) = xln 



where 

= _ M*) cos(x) 

For small x, i.e. i < 1, /(x) = 1/3 + 0(x 2 ). For typical physical parameters such that 
#| A | < 7r, we therefore observe a thermal-maser phase transition for values of 9 determined 
by Vq(x) = 0. The corresponding critical transition line is then 

a w=\ + \m^> (41) 

where ii(#,0) = Ii{9). Eq. ([Tj]) coincides with the radius of convergence of the thermal 
probability distribution Eq. ( |TTD as it should. For A = we see that 9 2 e ^ = + 9 2 and 
hence this phase transition can occur for 9 = 0, i.e. the phase transition can be induced 
by the presence of noise in the stochastic variable 9. For 9\/S\ ^> 7r we must in general 
make use of a combination of analytical and numerical methods, as in the case of ag = 
fT2|| , in order to get a detailed picture of the phase diagram. 

In the maser phase, the order parameter x(8) approaches zero when the system ap- 
proaches the critical line Eq. (f4"ip. Furthermore, in the large N limit x(8) is always zero in 
the thermal phase. Hence, the order parameter is continues on the critical transition line 
Eq. (f41"P. To determine the order of the phase transition we therefore have to investigate 
higher order derivatives. The first derivative of x(6) with respect to 9, x'(8), at the critical 
transition line Eq. ([y]) is 

x'(6) - (42) 
where I[{0) is the derivative of I\{9) with respect to 9 and where 

(43, 

By making use of the positive definiteness of P#(£)£ 2 , the Cauchy-Schwarz inequality 
implies |/2(#)| < V a e ^ ® 2 - ^ x'{9) is zero, then we have to investigate higher order 



1 + n b + a 9l ff 
n b + a9l ff 



+ nb(2a — 1) 
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derivatives. The second derivative of x(9) with respect to 9, x"(8), at the critical thermal- 
maser line Eq. ( f4~T|) is 

X[ "> h(9)-W) ' 1 ' 

where I['{9) is the second derivative of I\(9) with respect to 9 and i s the derivative 
of l2{9) with respect to 9. 

The effective potential Vq(x) approaches its large 9 limit exponentially fast as a func- 
tion of a 2 . This unique minimum is determined by dVo(x)/dx = 0, where dVo(x)/dx is 
given by Eq. (|33|) mthy(x) andw(x) replaced by y(x) = a e~ 2 ^ x+A ^ a e /[ § (l+e~ 2( - x+A2 '> CT o) + 
n b {x + A 2 ) ] and w{x) = b e -^+^ 2 W e /[ |(i + e -20r+A>f ) + ^ + n ^ x + A 2 ) ], respectively. 

If the micromaser is not detuned, i.e. A = 0, then there is one critical thermal- 
maser transition line only. This critical transition line can be found analytically for any 
probability distribution Pe{£,)- The integral I\{9) is then given by 

1^0) = a 2 + 9 2 >l . (45) 

The critical thermal-maser transition line is then only dependent of the variance a 2 and 
the mean value 9: 

^ = l +1 2^T¥- (46 > 

For values of a and 9 above this critical line, i.e. a > a(9), the micromaser system is in a 
maser phase. Below this line, on the other hand, the micromaser is in the thermal phase 
(see e.g. Fig. |). 

The order parameter x(9) is continuous on the critical line Eq. (|46|) . The first deriva- 
tive of x(9) with respect to 9, x'(9), at this critical thermal-maser line is 

^) = 6^y , (47) 

which is non-zero for Og < 1. The critical line Eq. (^) then describes a line of second-order 
(thermal-maser) phase transition. If a 2 > 1, Eq. (|46|) also describes a line of second-order 
phase transitions for all values of 9 except for 9 = where a(0) = 1/2 + l/(2a 2 ). Since 
x'{9) then is zero we have to investigate higher order derivatives. The second derivative 
of x(9) with respect to 9, x"(9), at the critical thermal-maser line is 
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Figure 2: The thermal-maser critical lines a{9) according to Eq. (^) for <r| = 
0.1,1,10, i.e. when the micromaser is not detuned. All other parameters are as 
in Fig. |]. When cig > 1 the critical thermal-maser line describes a third-order 
phase transition at a(0) = 1/2 + l/(2cr|) and a second-order transition for any other 
possible value of a. When 0% = 10, this third-order phase transition occurs at 
a(0) = 0.55. 



x"{9) 



(e 4 ) de + i2 



<e 6 > 



(4£ 



which is non-zero when the first derivative is zero. The point a(0) = 1/2 + l/(2erg) 
on the critical thermal-maser line therefore corresponds to a third-order transition. From 
Eq. (f48|) it follows, in addition, that we can have at most a third-order phase transition 
when fluctuations in 9 are taken into account. 

If the probability P#(£) is sufficiently peaked such that the its ^-variation is fast in 
comparison to ^-variation of sin 2 (£A), i.e. if <t | A| <C ir/2, then Eq. ([38D is reduced to 
I\{9) = sin 2 (0A ) /A 2 + cr 2 , cos (26* A ). The first thermal-maser critical line is then given 
by 
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1 1 A 2 
a{ ~ 9) = 2 + 2 sin 2 (#A) + a 2 A 2 cos(2#A) ' (49) 

which is consistent with the result of Ref. [|1^] for o 2 e = 0. 

If, on the other hand, the quantity P$(£) has a ^-variation which is slow in comparison 
to the ^-variation of sin 2 (f A) i.e. if cr e \A\ > vr/2, then Eq. fl3|) reduces to li(0) = 1/2A 2 . 



Eq. fl4T|) is then reduced to 

a(6) = 1 - + A 2 . (50) 

In order to get explicit analytic results we now choose the following gamma probability 
distribution for the pump parameter 

p o(0 = t^-ttC e-^ , (51) 
1 {a + 1) 

where f3 = 8/aj and a = 6 2 /a 2 — 1, so that (£) = 6 and ((£ — 6 1 ) 2 ) = a 2 . Other choices 
are possible, but are not, as long as the distribution is sufficiently peaked, expected to 
change the overall qualitative picture. The integral Ii(6,x) is then given by 



1 1 

2 x + A 2 



1 + 



2(z + A> 2 



cos I 77 arctan 



'2{x + A 2 )a 2 



e 2 /2 



<jz 



• (52) 



This averaged form of q(x) depends on the two independent variables 8 2 /(2<Jq) and 2(x + 
A 2 ) Oq only. In particular, if 



2ai 



(53) 



then Eq. (|52] ) reduces to 



/i(M) 



1 1 



1 - e- 2 ^ +A cos 29V^+A 



This explicit form of I\(9, x) can be achieved for a broad class of probability distributions 
e.g. a Gaussian distribution. If, in addition, 



1 



2a- 1 



(55) 
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Figure 3: The probabilities V{+) (solid line) and V(+,+) (dotted line) as a 
function of 6 when the q m -terms in Eq. (|2|) are averaged with respect to 6 with 
Ga = 25. All other parameters are as in Fig. [fl. In the large limit the probability 
V{+) approaches 7 ? o(+) = 0.5 and the probability V(+, +) approaches Voo(+, +) ~ 
0.38, which are indicated by solid lines in the figure. 



then any solution to the corresponding saddle-point equation 1 = (2a — l)Ii{0,x{6)) is 
exponentially close to 

x = a - 1/2 - A 2 . (56) 

If, on the other hand, 

^>0><r*>l , (57) 

then the corresponding saddle-point equation has one non-trivial solution only, i.e. Eq. (|56|) 
Under the conditions given by Eqs. (|53|) and (|55|) or Eq. (|57D the maser-maser phase 
transitions will be less significant. The phase diagram then essentially consists of critical 
thermal-maser line(s) only (see e.g. Fig. |^). With regard to the observable V(+), the 
requirement 9 > Q v= \ = 2t{N^/x + A 2 is compatible with Eq. (^3|) provided that iV is 
sufficiently large, i.e. 



16 




■ 
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6 

Figure 4: Comparison of the exact V(+) (solid line) and V{+) as given in 
Eq. (|24| ) (dashed line) with p(x) as given in Eq. ( |l8|) with x = 1/2. The parameters 
are the same as in Fig. ||. In this figure we only see the first two revivals, i.e. only 
the terms v = 0, 1, 2 in Eq. ( p4| ) is visible in this figure. 



iV»#^L, (5* 
2vr - 1/2 



for the saddle-point solution Eq. (p6|). The equilibrium probability distribution is in this 
case given by Eq. (|1^) with the variance a 2 x = [a + nb(2a — 1) ]/2N and Xj = x as in 
Eq. fl56|). The equilibrium distribution p(x) is peaked around x provided x ^ a x , i.e. 

1 q + n b (2a-l) 
iV 2 [a-l/2-A^p • [bJ) 

From Eq. (|27|) we observe that the u:th revival in V(+) occurs in the region where 9 is 
close to 



6 V = 2irvN^a - 1/2 . (60) 
According to Eq. (|2UJ), two consecutive revivals in V(+) are separated provided 
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v < V2N 



2a -1 1 



(61) 



y/a + n,b(2a — 1) 2 



which gives a bound on how many revivals that can be resolved in the excitation proba- 
bility V(+). In Fig. |3] we consider values of relevant physical parameters such that Eqs. 
(|55|), (|58|) and (|59| ) are satisfied. We then observe that V(+) and V(+,+) exhibit well 
pronounced revivals. For completeness we show in Fig. |], for the same set of values of 
the physical parameters as in Fig. ^ that the approximative but analytical expression 
for V(+) as given by Eq. (|24]) is in very good agreement with the exact result for V(+). 
Within one standard deviation we expect, according to Eq. flBTD , that at most seven re- 
vivals in V(+) are well separated in this case. Numerically we find that actually only the 
first two revivals are well separated. 

Let us now consider the case a\ <C 1. If, in addition, Eq. ( |53|) is satisfied, the quantity 
(q(x))g is reduced to Eq. ([I]). If on the other hand Eq. (|53|) is not satisfied, i.e. 6 <C of, 
then the saddle-point equation Eq. (P0|) has only the trivial solution for any probability 
distribution P#(£) as it should. The fluctuations in the atom cavity time have then no 
dramatic effect on the order parameter x (see e.g. Fig. |7]), the excitation probabilities 
V(+) and V(+, +) or the phase diagram. 



In this Section we consider the situation when there is a spread in the detuning param- 
eter A. Suppose that the spread in the detuning A is expressed in terms of a stochastic 
variable £ as described by the probability distribution Pa(0 such that (£) = A and 
((£ — A) 2 ) = a\. The averaged value of q(x) with respect to -Pa(£) is then 



4.2 Detuning Fluctuations 




(62) 



Non-trivial saddle-points of Vo(x) can then be found by solving the equation 



1 = (2a-l)Ji(0,x(0)) , 



(63) 



where 




(64) 
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For small values of x we follow the argument of Section f4.1| concerning the critical 
thermal-maser phase transition line. The critical transition line is then 



'^'l + l-W)- 1 ' (65) 



where Ji(#, 0) = Ji{6). Eq. fl6"5|) coincides with the radius of convergence of the thermal 
probability distribution Eq. (|TTD. In general we must use a combination of analytical and 



numerical methods, as in the case of ag = |[12|| , in order to get a detailed picture of the 
phase diagram. 

If the probability Pa{0 is sufficiently peaked such that the its £- variation is fast in 
comparison to ^-variation of sin 2 (#£), i.e. if a&6 <C ir/2, then Eq. fl6"4|) is reduced to 
Ji(0) = [3 sin 2 (#A)/A 2 + 6 2 cos(2#A) - 29 sin(2#A)/A ]/A 2 . The first thermal-maser 
critical line is then given by 

1 1 A 2 

<*(*) = » + 9 ■ 277TXTT 2 TflS . (66) 

2 2 sin [OA) + a A g{0) 



where g(9) = 3 sin 2 (6A)/A 2 + 9 2 cos(26A) - 29 sin(20A)/A, which is consistent with |L2 
for a\ = 0. 

If, on the other hand, the quantity -Pa(0 has a ^-variation which is slow in comparison 
to the ^-variation of sin 2 (#£), i.e. if a&6 3> ir/2, then J\[9) reduces to 



Ji(0) = 1/2 A 2 . (67) 



Eq. (|66"D is then reduced to 



a(#) = i + A 2 . (68) 

In order to get explicit analytic results we now choose the following exponential prob- 
ability distribution for the detuning parameter 

such that (£) = A and ((£ — A) 2 ) = cr A . Other choices are possible, but are not, as long as 
the distribution is sufficiently peaked, expected to change the overall qualitative picture. 
The averaged value of q(x) with respect to Pa(0 is given by 

(?(*))a = (q{x)) + (q{x)) osc , (70) 
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where 



It e a 

« w> ° s 27m 1* ^TF • (71) 



and 



( g (o;)) osc = / ^^^ cos ( 2 ^x + e 2 ) • (72) 



We observe that 



(9(^)0 = J £2 e-*l*Z+*l** I drj exp 



DC 



_ 2 1 A 2 X 
~ V + ?2oi2^i 



(73) 



Eq. ( |73|) can be used to find an expansion in the parameter A 2 /2a 2 x . For our purposes 
we notice that if o\ is sufficiently large, i.e. 

4»1 > (74) 



then Eq. (|73|) is reduced to (q(x)) < | a/tt x/2a\ + O(a;/2o" 2 1 ). Any solution of the 
corresponding saddle-point equation is much less then unity since (q(x)) osc < (q(x)) . 
The maser-maser phase transitions will therefore be less significant. The phase diagram 
does then essentially consist of critical thermal-maser line(s) only. For sufficiently large 
values of 6, i.e. 



the oscillating term can be estimated by \(q{x)) osc \ « [vr/8a|] 1/4 y/(2a - l)/(8a|#) for 
the solution 

^ =(2a -l) 2 -^ , (76) 

of the saddle-point equation. Eq. ( ff6|) is much less then unity due to the assumption 
Eq. (|74|). The term (q(x)) osc is negligible in comparison to (q(x))o when Eq. ([75]) is 
satisfied. With regard to the observable P(+), the requirement 9 > 6^=1 = 2tcN\ / x + A 2 
is compatible with Eq. (|75|) provided that iV is sufficiently large, i.e. 

N » ± sTt\-2 > ( 77 ) 
7r (2a — I)" 1 
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for the solution Eq. ([76]) of the saddle-point equation. The equilibrium probability dis- 
tribution is given by Eq. ( |I8D with Xj = x as in Eq. ( |76| ) and er 2 = (2a — l)[a + rib(2a — 
1) ]n/(4:Na\). The probability distribution p(x) is therefore peaked around x provided 
that x ^> a x , i.e. 

JV » 16|a + " t ' 2 °r 1)l «i • (78) 
(2a — !)•% 

The z/:th revival of P(+) occurs in the region where 9 is close to 

9 u = (2a-l)rrJ- — , (79) 

V 2 cr A 

according to Eq. (f2~7|). Eq. ([2^) implies that two consecutive revivals are separated pro- 
vided that 



N (2a -1)% 1 
U - X 'J[ 16[a + n 6 (2a-l)] ~2 ' (80) 

which gives a bound on how many revivals that can be resolved in the excitation proba- 
bility V(+). Fluctuations in the atom cavity transit time with a large aj give, however, 
an order parameter of the order one. This difference between 9- and A-ffuctuations is 
illustrated in Fig. [7J. In Fig. |5] we consider physical parameters such that Eqs. flT4]), ( |77| ) 
and ( |78"D are satisfied. We observe that V(+) and V(+, +) now exhibit well pronounced 
revivals. With the same set of physical parameters one can also verify that the approxima- 
tive but analytical expression for V(+) as obtained by a Poisson resummation technique, 
i.e. Eq. (|24|), is in good agreement with the exact expression for V(+). We also notice 
that according to Eq. (|80|) , only the first revival in V(+) is clearly resolved in this case, 
which agrees well with our numerical results. 

Let us now, on the other hand, consider the case when a\ is small but non-zero, i.e. 

< a 2 A < 1 . (81) 

The quantity (q(x))o can then be approximated according to (q(x))o = x/[2(x + A 2 ) ] and 
the term (q(x)) osc has an upper bound |(g(a;)) osc | < 1/[1 -l-e 2 ] 1 / 4 , where e = 2o\9x/{ x + 
A 2 ) 3 / 2 . For sufficiently large values of 9, i.e. 

1 (a-1/2) 3 / 2 , . 
g>> 2^ a-l/2-A 2 ' ^ 

the term {q(x)) osc is negligible in comparison to (q(x))o for the solution 
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Figure 5: The probabilities V{+) (solid line) and V(+,+) (dotted line) as a 
function of 9 when the g m -terms in Eq. (|2|) are averaged with respect to A with 
<7^ = 25 and N = 800. All other parameters are as in Fig. In the large 9 
and N limit the probability V{+) approaches T ) 00 {+) = 0.5 and the probability 
V(+, +) approaches 7 ? o(+ ; +) ~ 0.38, as indicated by a solid line in the figure. For 
sufficiently large 9, x w 0.02 and (V{+)) » « 1. 



x = a - 1/2 - A 2 , (83) 

of the saddle-point equation. The maser-maser phase transitions will therefore be less 
significant if Eqs. ([81]) and ([52]) are satisfied. The phase diagram does then essentially 
consist of critical thermal-maser line(s) only. For sufficiently small values of 9, i.e. 

t « ^ , (84) 

Eq. (|70l) is reduced to Eq. (|l|). Hence there is no dramatic change in the order parameter 
or in the phase diagram. With regard to the observable V(+), the requirement 9 > 9 V =\ — 
2nN\/x + A 2 is compatible with Eq. ( |5^ ) provided that N is sufficiently large, i.e. 
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Figure 6: The probabilities V{+) (solid line) and V(+,+) (dotted line) as a 
function of 9 when the g m -terms in Eq. (|2|) are averaged with respect to A with 
<7^ = 0.1 All other parameters are as in Fig. [fl. In the large 6 and N limit the 
probability V(+) approaches Voo(+) = 0.5 and the probability V(+, +) approaches 
7'oo(+)+) ~ 0.38, which is indicated by a solid line in the figure. For sufficiently 
large 8,iw 0.5, (V{+)) as 0.5 and {V {+,+)) ~ 0.38. 



1 Ja- 1/2 , . 

jV>> 4^i a-l/2-^ ■ (85) 

The equilibrium probability distribution is then given by Eq. flTB|) with Xj — X dbS 111 

Eq. ( |8"3"D and = [a + rib(2a — 1) ]/ (2N) according to Eq. (ITS). The distribution is 
peaked around x provided that x ^> a x , i.e. 

1 a + rih(2a — 1) , 

JV>> 2 a -l/2-^ ' (86) 
The w.th revival of V(+) occurs in the region where 9 is close to 



9 V = 2irvN Va - 1/2 , (87) 
according to Eq. (|27|) . Furthermore, two consecutive revivals are separated provided that 
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V < 



'N (2a -l) 2 1 



2 a + n b (2a-l) 2 ' 

which gives a bound on how many revivals that can be resolved in the excitation proba- 
bility V(+). 

In Fig. H we consider physical parameters such that Eqs. (|85| ) and (|86f) are 

satisfied. In this case we also observe well pronounced revivals in V(+) and V(+, +). The 
width and location of the revival in V(+) agrees well with the approximative expression 
Eq. (H. 



5 The Correlation Length 



Let us now consider long-time correlations in the large N limit as was first introduced 
in Ref.||. These correlations are most conveniently obtained by making use of the con- 
tinuous time formulation of the micromaser system fL0| . The lowest eigenvalue Aq = of 



the matrix L, as defined in Section then determines the stationary equilibrium solution 
vector p = p as given by Eq. (Q). The next non-zero eigenvalue A of L, which we deter- 
mine numerically, will then determine typical scales for the approach to the stationary 
situation. The joint probability for observing two atoms, with a time-delay t between 
them, can now be used in order to define a correlation length ■y A (t) |J. At large times 
t —>■ oo, we define the atomic beam correlation length £4 by || 

7A (t)~e-*/^ , (89) 

which then is determined by A, i.e. 7^4 = 1/A. For photons, we define a similar correlation 
length £c- It follows that the correlation lengths are identical, i.e. £a = £c* = £ 0- 

The correlation length 7^ is shown in Fig. [| for different values of <jg and exhibits large 
peaks for the pump parameters 0* kk+l , 9q and/or 9* k Ul2| . The correlation grows at most 
as \fN at 6 = 6q. At 6^ fc+1 and Q\ k the large iV dependence is, however, different. At these 
values of the pump parameter there is instead a competition between two neighbouring 
minima of V (x). Using the technique of 0] it can be shown that the peak close to 
= Otk+i is 
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Figure 7: The normalized steady-state mean photon number x = n/N as a 
function of 9 for different values of and a\ when a = 1, n& = 0.15, A = and 
iV = 1000. In the large and iV limit the order parameter approaches a constant 
value. This asymptotic value is determined by Eq. (33) when a 2 e = o\ = 0. The 
solid line at ps 0.34 in the figure corresponds to this asymptotic value of the 
order parameter. 
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7 £^ e ^> , (90) 

where AVo is the smallest potential barrier between the two competing minima of the 
effective potential V (x). For increasing values of aj or a\ this potential barrier AVo 
decreases. The corresponding correlation length will then also decrease (see e.g. Fig. H). 
Numerical study shows that adding fluctuations to the atom-photon frequency detuning 
result in the same quantitative effect: the peaks in the correlation length decreases due 
to smaller potential barriers AVq. 



6 Conclusion 
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In conclusion, we have investigated various effects of including noise-producing mech- 
anisms in a micromaser system, like a velocity spread in the atomic beam or a spread in 
the atom-photon frequency detuning. For sufficiently large widths of the corresponding 
probability distributions, i.e. for sufficiently large a% or a non-zero a\, the excitation 
probabilities of atoms leaving the microcavity exhibit well-pronounced revivals. Noise 
therefore tends to increase the order in the system. Furthermore, the maser-maser phase 
transition disappear for sufficiently large aj or a\, in which case the phase diagram in the 
a- and 9- parameter space consists essentially of a single maser-thermal critical line. We 
have also shown that the correlation length drastically decreases when noise-producing 
mechanisms are included in the system. 
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APPENDIX 



In this Appendix we derive an exact analytical expression for Vq(x) in the large 9 limit 
when the quantity q n has no spread in any of the physical parameters, i.e. when q(x) is 
as given in Eq. (P. The potential Vq(x) may then be re- written in the from 



V {x) 



dul In 



v 



+ A 2 



+ n b v 



+ In 



1 + 



,sin 2 (# v / ^ + A 2 ) - 
a + n b (u + A 2 ) 



(91) 



In 



v + A 2 



+ {l + n b )v 



In 



1 + 



b sin 2 (0\/^ + A 2 ) - b 
b+ (1 +n b )(v + A 2 ) 



By making use of the power series expansion of the logarithm we obtain 



VJx) 



- J dv{ hi 



& — — — + (1 + n b )v ) / ( a — + n b v 



v + A 2 



z/ + A 2 



(92) 



+ E 



^ nL V.a + n 6 (z/ + A 2 )/ V & + (1 + 00 + A2 ) 



cos 2n (^vVTA 2 ) 



In the large 6 1 limit only the non-oscillating terms in Eq. fl92|) contribute, i.e. 







"1 - a 




In 






a 



+ f(y(„))-f( W (y)) 



(93) 



where 



f{z)=\nz + R{z) , (94) 

and y(x) = a/(a + n b (x + A 2 ) ) and w(x) — b/( b + (1 + n b )(x + A 2 ) ). The function R(z) 
is defined by 



^ 2^ 2 2 ™(n!) 2 ' 

n=l v 7 



(95) 



which is a generalized hypergeometric series JTj]. By differentiating Eq. (|93|) we therefore, 
finally, obtain Eq. (|33|). 



28 



ACKNOWLEDGEMENT 



The authors wish to thank the members of NORDITA for the warm hospitality while 
the present work was completed. The research has been supported in part by the Research 
Council of Norway under contract no. 118948/410. 

References 

[1] Y. Yamamoto and H.A. Haus "Preparation, Measurement and Information Capacity 
of Optical Quantum States", Rev. Mod. Phys. 58 (1986) 1001; Y.I. Vorontsov "Stan- 
dard Quantum Limits of Measurement Errors and Methods of Overcoming Them" , 
Sov. Phys. Usp. 37 (1994) 81; CM. Caves and P.D. Drummond "Quantum Limits 
on Bosonic Communication Rates" , Rev. Mod. Phys. 66 (1994) 482; C.H. Henry and 
R.F. Kazarinov "Quantum Noise in Photonics", Rev. Mod. Phys. 68 (1996) 801. 

[2] A. Maritan and J.R. Banavar, "Chaos, Noise and Synchronization" , Phys. Rev. Lett. 
72 (1994) 1451; L. Gammaitoni, P. Hanggi, P. Jung and F. Marchesoni "Stochastic 
Resonance", Rev. Mod. Phys. 70 (1998) 223; J.W.C. Robinson, D.E. Asraf, A.R. 
Bulsara and M.E. Inchiosa, "Information-Theoretic Distance Measures and a Gener- 
alization of Stochastic Resonance", Phys. Rev. Lett. 81 (1998) 2850; H. Zhonghuai, 
Y. Lingfa, X. Zuo and X. Houwen, "Noise Induced Pattern Transition and Spa- 
tiotemporal Stochastic Resonance" , Phys. Rev. Lett. 81 (1998) 2854; G. Giacomelli, 
F. Marin and I. Rabbiosi, "Stochastic and Bona Fide Resonance: An Experimental 
Investigation", Phys. Rev. Lett. 82 (1999) 675; S. Kim, S. Hee Park and H.-B. Pyo, 
"Stochastic Resonace in Coupled OScillator Systems with Time Delay 1 , Phys. Rev. 
Lett. 82 (1999) 1620; D. Nozaki, D.J. Mar, P. Grigg and J.J. Collins, "Effects of 
Colored Noise on Stochastic Resonance in Sensory Neurons" , Phys. Rev. Lett. 82 
(1999) 2402; T. Ohira and Y. Sato, "Resonace with Noise and Delaf , Phys. Rev. 
Lett. 82 (1999) 2811; V.S. Anishchenko, A.B. Neiman, F. Moss and L. Schimansky- 
Geier "Stochastic Resonace: Noise- Enhanced Order 11 , Sov. Phys. Usp. 42 (1999) 7; 
Y. L. Klimontovich "What are Stochastic Filtering and Stochastic Resonance? 1 , Sov. 
Phys. Usp. 42 (1999) 37. 

[3] A. Buchleitner and R.N. Mantegna, "Quantum Stochastic Resonance in a Micro- 
maser 11 , Phys. Rev. Lett. 80 (1998) 3932. 



29 



[4] H. Walther, "The Single Atom Maser and the Quantum Electrodynamics in a Cav- 
ity" , Physica Scripta T23 (1988) 165; "Experiments on Cavity Quantum Electrody- 
namics" , Phys. Rep. 219 (1992) 263; "Experiments With Single Atoms in Cavities 
and Traps" in "Fundamental Problems in Quantum Theory 1 , Eds. D.M. Greenberger 
and A. Zeilinger, Ann. N.Y. Acad. Sci. 755 (1995) 133; "Single Atom Experiments 
in Cavities and Traps" , Proc. Roy. Soc. A454 (1998) 431; "Quantum Optics of a 
Single Atom", Laser Physics 8 (1998) 1 and in Physica Scripta T76 (1998) 138. 

[5] B. W. S. and P. L. Knight, "The Jaynes-Cummings Model", J. of Mod. Opt. 40 
(1993) 1195. 

[6] "Nonclassical Effects in Quantum Optics", Eds. P. Meystre and D.F. Walls (AIP, 
1991). 

[7] B.-S. Skagerstam in "Applied Field Theory" , Eds. Choonkye Lee, Hyunsoo Min and 
Q-Han Park (Chungbum Publ. Hose, Seoul, 1999). 

[8] D. Filipowicz, J. Javanainen and P. Meystre, "The Microscopic Maser", Opt. Comm. 
58 (1986) 327; "Theory of a Microscopic Maser" Phys. Rev. A34 (1986) 3077. 

[9] P. Elmfors, B. Lautrup and B.-S. Skagerstam, "Correlations as a Handle on the 
Quantum State of the Micromaser" , CERN/TH 95-154 ( |cond-mat/9506058|) ; Physica 
Scripta 55 (1997) 724; "Correlations in the Micromaser" Phys. Rev. A54 (1996) 
5171. 

[10] L. Lugiato, M. Scully and H. Walther, "Connection Between Microscopic and Macro- 
scopic Maser Theorf, Phys. Rev. A36 (1987) 740. 

[11] E.T. Jaynes and FW. Cummings, "Comparison of Quantum and Semiclassical Ra- 
diation Theories With Application of the Beam Maser" 1 , Proc. IEEE 51 (1963) 89. 

[12] P.K. Rekdal and B.-S. Skagerstam, "On the Phase Structure of the Micromaser 
System", Opt. Commun. 184 (2000) 195 and "Theory of the Microscopic Maser 
Phase Transitions" , CERN-TH/99-410 (submitted for publication). 

[13] M. Fleischhauser and W.P. Schleich, "Revivals made simple: Poisson summation for- 
mula as a key to the revivals in the Jaynes-Cummings model" , Phys. Rev. A47 (1993) 
4258; R. Courant and D. Hilbert, "Methods of Mathematical Physics 11 , Intersience, 
New York, 1953. 



30 



A. Erdelyi, W. Magnus, F. Oberhetter and F. G. Tricomi, "Higher Transcendental 
Functions" , McGraw-Hill Book Company, Volume I, 1953. 



31 



